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ABSTRAGT 

We present a specific periodicity of the Galaxy images in the multipole domain which can be used 
for separation of the cosmic microwave background (CMB) signal from the Galaxy and foreground 
component. This method takes into consideration all the coefficients of the expansion of the signal 
from the sky AT(0,0) into spherical harmonics and dealing with combination of multipoles ai^rn and 
<^z+A,m for each multipole mode (^, —i<m<i) from the whole sky without galactic cut, masks or 
any dissection of the whole sky into disjoint regions. For the polar coordinate system we use particular 
values for A = 4n, n = 1, 2,... which remove all bright point-like sources localized in the Galactic plane 
and strong diffused component down to the GMB level. To illustrate that significant correlations of the 
point like sources and foregrounds coefficients we apply this method to the WMAP Q,V 

and W bands in order to remove these marks of Galaxy from the maps. We believe that our method 
would be useful for separation of the CMB signal from different kind of ’’noises” in the maps. 


1. INTRODUCTION 

Separation of the CMB signal and foregrounds (espe¬ 
cially, Galactic foregrounds) is the major problem for 
all the CMB experiments, including ongoing the WMAP 
and future the PLANCK missions. The basics approach 
to solve this problem is related with implementation of 
different masks and disjoint regions of the whole sky 
maps (see,Bennett et al. (2003a),Bennett et al. (2003b) 
Tegmark, de Oliveira-Costa & Hamilton (2003),Eriksen et 
al. (2004) ). However, the question is how the Galactic re¬ 
gion of the map can affected any estimators of the primor¬ 
dial CMB signal, particularly, the power spectrum C{1) 
and phases of the multipole representation of the CMB, if 
we try to process the data sets without any masks? This 
question closely related to the question, what statistical 
properties of the Galactic foreground could be? 

In this letter we would like to present some of the results 
of investigation of the Galactic area of the WMAP sky us¬ 
ing the K — W bands signals. Particularly, we will show 
that in the polar coordinate system the brightest area of 
the map related with the Galaxy can be removed down 
to the CMB level by simplest linear combination of the 
spherical harmonics coefficients ai^m 

lA _ /-.x 

^Lm — ^l,m I I RZ+A,m Gj 

|R/+A,m I 

where A is free parameter, which corresponds to transi¬ 
tion from the I multipoles to / + A multipoles for the same 
\m\ < I- modes. As one can see from Eq. (1) our idea is to 
introduce new characteristics of the AT(^,0) field on the 
sphere 

^max I 

= E E dtmYi,m{e,cp) (2) 
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(here Yi^rniO^ (j) are spherical harmonics) which can signifi¬ 
cantly decrease the Galactic plane signal in the K-W maps. 
In Section 1 we show, that this new estimator directly 
related with phases and power spectrum of the whole sky 
signal and can be used for separation of the CMB signal 
from foregrounds. Section 2 is devoted to practical imple¬ 
mentation of the method to the WMAP data sets. In this 
section we will show that for with A = 4n,n = 1,2.. 
the D^iO^cj)) map derived from K-W maps do not have 
strong signal from the Galaxy. In Section 3 we discuss ro¬ 
tational invariance of A = 4n correlations and show that 
they are specific for particular orientation of the Galaxy in 
the reference system. In the conclusion we summarize all 
the results and briefly discuss an application of the method 
for cleaning of the CMB maps from foregrounds. 

2. PROPERTIES OE ESTIMATOR. 

To understand the properties of d^^ estimator let us de¬ 
fine the transformation of the AT(0, 0) signal to spherical 
harmonic coefficients ai^rn RS 

Imax I 

AT(0,<^)=E E (3) 

1=2 m=—l 

where \ai^rn\ and ^i^rn are the magnitudes and phases of 
/, m harmonics correspondingly. In terms of these variables 
our estimator d^^ now has a form 

dtm = \ai,mT‘^‘'"' (1 - - i'i+A.m]) (4) 

One can see, that d^^ coefficients depend on the phase 
difference taking for /, / + A multipoles. If 

the phase difference is small enough ^i^rn — ^/+A, <C 7r/2, 
then 
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So, for Gaussian CMB signal the phase difference is not 
especially small because of non- correlation of the phases, 
while for non-Gaussian signal these correlation should be 
significant (see, for example, Chiang et al, 2003, Naselsky 
et al, 2004). According to Eq. (5) their contribution to 
d^m became to be vanish, if is almost the same as 

To illustrate that tendency, let us describe the simplest 
model (see Naselsky et al, 2004), when a bright point 
sourse is placed in the Galactic plane area, having coordi¬ 
nates Oj = 7 r/ 2 , 0 j and amplitude Aj. The phases of the 
harmonics for that signal are (Naselsky et al, 2004) 


tanT/,^ = - 


sin(m0j) 

cos(m0j) 


(6) 


and they do not depend on /. So, for single point source 
^i,m = ^z+A,m and that source do not contribute to the 
D^(0, (p) map for all A = 1 , 2 ,... in Eq(l). Eor combina¬ 
tion of the sources AT(^, 0) = ATj{0, (j)) the properties 
of the phases are more complex that for the single one (see 
Naselsky et al, 2004), but the tendency of ’’selfcleaning” 
of the D^{0^ (j)) map should be preserved, at list for some 
values of A-parameter. In the next section we will show 
that for A = 4n, where n = 1 , 2 .. this method provide the 
best result for the Galactic plane area, including not the 
point sources only, but diffused foreground as well. 


3. ”4n” periodicity of the wmap bands. 

To investigate the properties of estimator we use the 
HEALPix package (Gorski,Hivon & Wandelt, 1999) to de¬ 
compose each of the WMAP signal maps at K-W bands for 
the spherical harmonics coefficients (Wide = 512). 

Then we perform transformation of the coefficients 

to d^^ for each band and obtain the following maps, 
shown in Eig.l. We reproduce the same analysis using 
the GLESP package (Doroshkevich et all, 2003) in order 
to show that corresponding properties of the maps do not 
depend on the sky pixelization . 

To remove contribution of the GMB signal to the 
D^{0^(j)) map, in Eig.2 we plot for difference 

between V and W maps, which containt the foregrounds 
and instrumental noise only. 

In Eig.3 we show the D^iO^cj)) maps for W band and 
different Imax^ but for the same value A = 4. Note that 
for Imax > 150 the properties of 0) map depend 

on the instrumental noise, while for Imax = 50 all the 
2 )A= 4 (^, 0 ) have quite regular structure, without 

Galactic belt. Important remark is that if A 7 ^ 4n then 
d^^ estimator does not remove the strong signal in the 
Galactic plane. 

4. ROTATIONAL INVARIANCE OF THE d'^"^(6>,0) MAP . 

The estimators of Galactic signal are dealing with 
phases of ai^rn coefficients and consequently depend on the 
reference system of coordinates . Obviously, for differ¬ 
ent system of coordinates these ai^rn coefficients and cor¬ 
responding phases are different for the same values l^m. 
Practically speaking, are estimators rotationally in¬ 

variant and if not, how significantly that non-invariance 
can transforms any conclusions about contribution of the 


Galactic signal? To answered to these questions we need 
to know how to transform a given set of ai^rn taking 
for given reference coordinate system to new coefficients 
bi^rn which corresponds to new coordinate system rotated 
by Euler angles Eollow to general method Var- 

shalovich,Moskalev & Khersonskii (1988) and taking into 
account properties of spherical harmonics (see Goles ,2005, 
Naselsky et al, 2004) we get 

bl,m = E /5,7)az,„j' (7) 

m' 

where (a,/3, 7 ) is a spherical harmonic decom¬ 

position of the Wigner function D((T,/?, 7 ). The co¬ 
efficients (a,/3, 7 ) should preserve the modules 

under transformation. With¬ 
out lost of generality let assumes that for initial ref¬ 
erence system of coordinates ai^rn = exp{i^i^rn) where 
are the phases for given I and \m\ < 1. To ob¬ 
tain the trigonometric moments in the reference sys¬ 
tem after rotation we define a matrix of correlations 



Fig. 1.— From the top to the bottom : the (p) map for 

the V band for A = 0 , A = 1 and A = 4. Next tree maps are the 
same , but for W band. All the maps are taking for I < 250. Here 
and bellow the amplitudes of the signal are in micro Kelvin. For all 
^A= 4 (^ 0 ) jnaps tjie range of temperature scale is —0.5,0.hmAT. 
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where I = / + A, m = m and bi^rn = \bi,m\exp{i(l)i^rn), 
(j)i^rn) is the phase of l,m harmonics after rotation. So, 
under rotation the (p) map transforms by the fol¬ 

lowing way: 


-^m,p(«> /3, (3, 

p,q 


(9) 

Taking into account that (see Varshalovich,Moskalev & 
Khersonskii (1988)) 

^L,m’ = exp(-iTO'7); (10) 


where 


m' (/5) = {I + m )\{l - m )\{l + m)\{l - m)\ 


min(l—m,l—m ) 

E (-1)' 

k 


—m —k 


{I — m — k)\k\ 


{I — m' — k)\{k + m + m')\ 


X (cos (sin 

( 11 ) 


we obtain the following formula : 


where 


/ ^ m,m ,p,q ^ ' 


W^’^' , (/3) = „(/?)/, (/3) 


( 12 ) 


The Eq(12) show that after rotation the properties of 
the phase correlation became to be more complicated than 
in the reference system. To estimate their analytically, we 
would like to point out that real and imagenary parts of 
the exponent in Eq. (12) now containts not phase dif¬ 
ference — ^i' only, but additional terms j{p — q) 

and a{m — m). However, taking into account that in the 
reference system the major correlations corresponds to 
q = p and I = / + 4n, we can draw our attention on these 
modes only. As one can see from Eq. ( 12 ), for these modes 
dependency of exponent in the left hand side on a and 7 
parameters is vanished , while dependency on [3 looks like 
a linear combination of exponents — ^ij^au.p convolved 
by the coefficients Obviously , that part of 

exponent is not rotationally invariant, while rotation of 
the map with /3 = 0 preserved rotationally invariance of 
the phase difference. To illustrate that effect in Eig.4 we 
plot the T)^^^(^,0) map taking from the W band, and 
rotated by the angles 0 = 77"^, p = 38*^ (which corresponds 
to 7 ^ 0 ) relatively to the system of standard polar coor¬ 
dinates. As one can see from this fig. the Galactic belt is 
still the major source of the high amplitude signal. Thus, 
the effect of 4n-periodicity is rotationally invariant for par¬ 
ticular choosing of the Euler angles a ^ ^ 0, f3 = 0. 




Fig. 2.— The maps (p) for the difference V -W bands 

for A = 4 and Imax = 50(top). Bottom - the same map, but for 
Imax — 150 





Fig. 3.— From the top to the bottom : the (p) maps for 

W band for Imax = 150 and Imax = 50. 



Fig. 4.— The cp) map taking after rotation of the W 

band map by the angles 6 = 77°, (p = 38°. . 

At the end of this section we would like to point out 
that for given statistical properties of the primordial 
CMB signal the phases of estimator and power 

spectrum closely related to the phases ^i^rn and power 
spectrum of the signal C{1), obtained using airn co¬ 
efficients. Eor primordial Gaussian signal these rela¬ 
tionships are extremely simple. The phases of 
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estimator should be non-correlated and uniformly dis¬ 
tributed within interval [0,27r], while the power spectrum 
= 1 /( 2 / H- 1 ) related to standard one 

C{1) = 1/(21 + 1 ) Y,m through equation 

^ m=l 

D^=\l) = 2C{1) - ^ \ci,rnf COs{^l+A,m - £.l,m) 

m=—l 

From Eq. (13) one can see that ^ 2C{1) with 

error about ’’cosmic variance” error. (^ } ^ ). For non- 

V^+ 2 

Gaussian primordial CMB signal with correlated phases 
6 +A,m 56 ,m the power of signal has a tendency to 

drop down to the level 0 , if ^z+A,m ^ 6 ,m- 

That tendency open a new possibility to test Gaussian- 
ity of the GMB signal trough power spectrum estimation, 
using aim and coefficients and corresponding power 

spectra C{1) and For Gaussian signal we should 

obtain a/ 2(7(/), while for non-Gaussian signal it 

should be < 2C{1) with ’’cosmic variance” error 

bars. 


5. CONCLUSION. 

In this paper the new estimator for the GMB maps has 
been presented. This estimator is based on the phase 
correlation between / + 4n, m and /,m harmonics and it 
decrease significantly contribution of the Galactic fore¬ 
grounds to the map. Decrovered 4n, n = 1 , 2 ... periodicity 
of the Galactic foegrounds phases seems to be important 
for accurate solution of the problem of separation the pri¬ 
mordial GMB signal from the foregrounds and instrumen¬ 
tal noise including whole sky signals. This method will be 
publish in a separate paper. 
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